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We have calculated the temperature dependence of shear 77 and bulk £ viscosities of quark matter 
due to quark-meson fluctuations. The quark thermal width originating from quantum fluctuations 
of quark-7r and quark-cr loops at finite temperature is calculated with the formalism of real-time 
thermal field theory. Temperature-dependent constituent-quark and meson masses, and quark- 
meson couplings are obtained in the Nambu-Jona-Lasinio model. We found a non-trivial influence 
of the temperature-dependent masses and couplings on the Landau-cut structure of the quark self¬ 
energy. Our results for the ratios 77 /s and £/s, where s is the entropy density (also determined in the 
Nambu-Jona-Lasinio model in the quasi-particle approximation), are in fair agreement with results 
of the literature obtained from different models and techniques. In particular, our result for rj/s has 
a minimum very close to the conjectured AdS/CFT lower bound, 77 /s = l/An. 
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I. INTRODUCTION 

A high-temperature, weakly interacting medium is nat¬ 
urally expected to be produced in heavy ion collision 
(HIC) experiments at high energies because the con¬ 
stituents of the medium (quarks and gluons) should be¬ 
come almost free according to the asymptotic freedom 
property of quantum chromodynamics (QCD). However, 
the experimental finding of elliptic flow in HICs at the 
Relativistic Heavy Ion Collider (RHIC) leads to the inter¬ 
pretation that the medium produced in the collisions is 
actually a strongly interacting liquid, instead of a weakly 
interacting gas. This interpretation comes, in first place, 
from the ability of hydrodynamics to describe the RHIC 
data with a small value of the shear viscosity over en¬ 
tropy density ratio, 77 /s M- The shear viscosity coef¬ 
ficient 77 of a medium represents the ability of its con¬ 
stituents to transfer momentum over a distance compa¬ 
rable to their mean free path. When the inter-particle 
coupling is strong, momentum transfer takes place eas¬ 
ily and, therefore, the shear viscosity of the matter is 
small. Several theoretical studies have estimated the ra¬ 
tio 77 /s, at weak and strong couplings and for different 
temperature regimes — a set of such studies is repre¬ 
sented by Refs. [914361 [3814421 . For earlier studies, see 
e.g. Refs. [43 46j. An interesting feature of the results in 
Refs. [TBl fl 6 l. I 21 L 28i] is that 77 /s reaches a minimum in the 
vicinity of a phase transition; the smallness of this mini¬ 
mum is significant in connection with a conjectured lower 
bound, 77 /s = 1/47T, known as the KSS bound, obtained 
in the context of AdS/CFT correspondence |47j . 

Similar to the shear viscosity, another transport coef¬ 
ficient is the bulk viscosity, £. It is defined as the pro¬ 
portionality constant between the non-zero trace of the 
viscous stress tensor to the divergence of the fluid ve¬ 
locity, and usually it appears associated with processes 
accompanied by a change in fluid volume or density. The 


bulk viscosity has received much less attention than the 
shear viscosity in hydrodynamical simulations because its 
numerical value is assumed to be very small, as it is di¬ 
rectly proportional to the trace of the energy-momentum 
tensor which generally vanishes for conformally symmet¬ 
ric matter. However, lattice QCD simulations [48j, l49j 
have shown that the trace of the energy momentum ten¬ 
sor of hot and dense QCD might be large near the QCD 
phase transition, which in turn indicates a non-zero, and 
possibly large value of £ as well as of £/s near the tran¬ 
sition temperature. In the framework of pure-gauge lat¬ 
tice QCD, C/s near the transition temperature is esti¬ 
mated in Ref. [so| . Analytical calculations employing 
different techniques and models have been used for esti¬ 
mating C of strongly interacting matter, a list of refer¬ 
ences are [26|, [H, IdTITgI] ; some of those [5ll453j indicate a 
divergent behavior for C near the transition temperature. 

In this work we em ploy the two-flavor Nambu-Jona- 
Lasinio (NJL) model [ 63 ] to evaluate the temperature 
dependence of the shear and bulk viscosities in the vicin¬ 
ity of the crossover temperature. The NJL model is a 
very useful model for studying many aspects of the chi¬ 
ral structure of QCD in vacuum and at finite temperature 
and baryon density [63]. The model has also brought a 
great deal of insight into the problem of viscosities of 
strongly interacting matter [26|, Ijfl HT gL [j2l . l45l . l6ll |. A 
practical and transparent way to evaluate the contribu¬ 
tions of quark-meson fluctuations to viscosities is to eval¬ 
uate one-loop quark self-energy diagrams to obtain the 
quark relaxation time or, equivalently, the quark thermal 
width [H|36|. Adopting this method, we have first an¬ 
alyzed the detailed Landau cut structures of the quark 
self-energy for quark-7r and quark-cr loops, where the tem¬ 
perature dependence of the quark and meson masses and 
their couplings play an important role for its on-shell 
contributions. The present work explicitly demonstrates 
the nontrivial contributions of the temperature depen- 
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dence and momentum dependence of the quark thermal 
width to the viscosities, in the kinematic domains where 
the quark pole remains within the regions of the Lan¬ 
dau cut. The new aspect explored in the present work 
relates to the contributions of quark-er and quark -tt fluc¬ 
tuations to 77 and £, which become the main sources of 
dissipation beyond temperatures where pion decay into 
on-shell quark-antiquark pairs becomes possible. This 
temperature is commonly known as Mott temperature 
(7 m) above which the threshold condition of pion dis¬ 
sociation i.e. m 7r ( T) > 2 Mq{T) is satisfied. In a very 
recent work (64J, the NJL model was used to investigate 
the role played by quark-meson fluctuations on the shear 
viscosity. Our work is complementary to that one, as we 
use a different formalism to evaluate viscosities, consider 
also the bulk viscosity, use different quark- 7 r and quark-er 
couplings reflecting dynamical chiral symmetry breaking, 
and make a detailed study of the Landau cut structure of 
the quark self-energy. The paper is organized as follows: 
in the next section, we address the formalism of com¬ 
puting shear and bulk viscosities in terms of the quark 
thermal width, which is deducted from the quark self¬ 
energy diagram in the framework of real-time thermal 
field theory (RTF). Next, in results section, the analytic 
structure of the quark self-energy is rigorously discussed 
and its contribution to shear and bulk viscosity coeffi¬ 
cients of quark matter is also discussed. A summary and 
conclusion is presented in the last section. 


II. FORMALISM 


In the Kubo formalism [65|, [ 66 j, the shear and bulk 
viscosities are defined in the Lehmann spectral represen¬ 
tation of the two-point correlation functions of operators 
involving the components of the energy-momentum ten¬ 
sor [67| : 

(")= lie, (1) 

\Cj *.l«il-o+ q° \ \A c {q\q) ) 


with the spectral functions A v (q) and A^(g) given by 


A v (g) = 

1 d 4 ®e < «- B ([7r«(s),7r«(0)]) j9 , 

(2) 

Ac(g) = 

1 d 4 xe^([P(x),V( O)])* 

(3) 


where 

Tr^Or) = T v(x)-^T k k , (4) 

V{x) = -br\(x)~clT m {x) . (5) 


description of the medium, an approximation we use in 
the present paper — a clear discussion on this issue can 
be found in Ref. 0 - 

We use the NJL model to obtain the above correlation 
functions. The Lagrangian density of the model for u 
and d flavors is of the form 

C = - m Q )ip + G [(4>4>) 2 + (tpi^Tip) 2 ] , (6) 

where mg = ( m u ,md ) is the current quark mass ma¬ 
trix and r = (r 1 ^ 2 ^ 3 ) are the flavor Pauli matrices. 
The energy-momentum tensor is given in terms of the 
Lagrangian density by 

dC — 

T^ = -^£ + __= + # (7) 

o{a^ip) 

The model will be treated in the quasi-particle approxi¬ 
mation or, equivalently, in the leading-order approxima¬ 
tion in the 1 /N c expansion, where N c = 3 is the num¬ 
ber of colors, which is also equivalent to the traditional 
Hartree approximation of many-body theory. In this ap¬ 
proximation, the thermal spectral functions can be ex¬ 
pressed in terms of the quark propagator. 

We use the formalism of real-time thermal field theory 
(RTF) to obtain the correlation functions. In RFT, the 
two point function of any field-theoretic operator has a 
2 x 2 matrix structure reflecting the time ordering with re¬ 
spect to a contour in the complex plane [fill. The relevant 
matrix can be diagonalized in terms of a single analytic 
function, which determines completely the dynamics of 
the corresponding two-point function. In particular, the 
retarded correlation functions needed for the evaluation 
of A v (q) and A^ (q) can be written in terms of of the 11 
component of the corresponding two-point functions ■ 
see Ref. [69j] for details. For example, ignoring for the 
moment quark-meson fluctuations, A v (q) can be written 
as 


A v (q) = 2 tanh 



Im Ilii(g), 


( 8 ) 


with 


U n (q)=iN c N f J ^ N(q,k)D 11 (k)D 11 (q + k), (9) 

where Dn{q) is the scalar part of the 11 component 
quark-propagator matrix (in the zero width case): 


Dn(fc) = 


-1 


k o ~ H ) 2 + ie 

—2-iriujQ n Q {k)8{kl - (w^) 2 ) . (10) 


In the above equations, c 2 is the sound velocity and ((..))p 
denotes an appropriate thermal average. We note here 
that the second term in the expression for V(x) is neces¬ 
sary to account for energy conservation in a quasi-particle 


The nQ{k) in Eq. IIUI) is the Fermi-Dirac distribution 
n Q (k) = 1 


1 + ’ 


( 11 ) 
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with u k Q = (fc 2 + Mg) 1 / 2 , and 


N{q,k) 


32 


3 

-4 


ko{k 0 + qo)k ■ (k 


k ■ (k + q) + ifc 2 (fc + q) 2 . (12) 
o 


Fig. DJa) shows this quark-quark loop diagram, which 
can be considered as a schematic representation of shear 
viscosity coefficient at the zero frequency and momentum 
limit. The quark mass Mq is the solution of the gap 
equation 


for the 7r coupling, and 

£QQa = gQQaipv^, ( 17 ) 

for the a coupling. Note that the couplings gQQ a and 
SQQtt are not bare couplings; they incorporate the effects 
of dynamical chiral symmetry breaking and as such are 
different from each other. At high temperatures, when 
chiral symmetry is restored, they become equal to each 
other, as we discuss in the next section. The meson 
masses tom and quark-meson couplings gQQM are calcu¬ 
lated from well-known expressions of the NJL model : 

1 — 2GIIm(w 2 = m 2 M ) = 0, (18) 


Mq = rriQ + 4N c NfGMQ 


d 3 k 1-2 n Q {k) 

(2tt)3 ■ ( } 


gQQM — 


dn M (co 2 ) 


duj 2 


(19) 


The temperature-independent part of the integral above 
is ultraviolet divergent, while the temperature-dependent 
part, which contains the Fermi-Dirac distribution func¬ 
tion n,Q(k), is finite. We use a ultraviolet cutoff A to 
regularize the divergent integral, with A fitted (together 
with the other parameters of the model, toq and G) to 
obtain reasonable values for the quark condensate, pion 
leptonic decay constant and the pion mass in vacuum. 
The finite integral is integrated without a cutoff. 

To proceed with the evaluation of the viscosities, we 
include the dissipative processes due to quark-meson fluc¬ 
tuations. The fluctuations introduce an imaginary part 
in the quark self energy giving a thermal width Tq to 
the quarks; as it stands, with no finite imaginary part in 
the quark propagator, Eq. ([9| leads to divergent viscosi¬ 
ties. Adding an on-shell, momentum-dependent thermal 
width TQ^k) in the quark propagator, one obtains [6f| 
for the shear viscosity the expression 

= 4 N c N f f d 3 k k A n Q (k) [1 - n Q (k)} 

V 15 T J (2tt) 3 (^) 2 r 0 (fc) ' 1 j 

This expression for rj is identical to the one obtained from 
the standard relaxation time approximation [28l . 143} . For 
bulk viscosity, from Eq. m we obtain: 


n M (w 2 ) = 2N c N f f d ' k l 2n k Q{k) F m (u 2 ), (20) 

J (Z7TJ UJq 


with 


TMw 2 ) = 
F«{u 2 ) = 


H) 2 

( w q) 2 — w 2 /4’ 

(21) 

^ IT 

OJOs- 

1 

(22) 


where Mq is the solution of the gap equation in Eq. (fl3l) . 
The integrals in Eq. (TSUI) are to be understood as 
principal-value integrals when w 2 > 4 Mq. Like with the 
gap equation, the temperature independent part of the 
integral for II m{u 2 ) is divergent, and we use the same 
cutoff A as in the gap equation. 


Q (k) 



(M=n, a) 
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4 N c N f 
T 



d 3 k n Q {k) [1 - n Q (k)] 


(2tt)3 


K) 2 r Q (fc) 



(15) 


which is also identical to the one obtained in Ref. [28} 
which uses the relaxation time approximation. 

We evaluate the quark thermal width TQ(fc) from the 
quark-meson ( QM ) loop contributions to the quark self- 
energy, pictorially represented in FigOJb). We need the 
meson masses and quark-meson couplings. We assume 
the couplings of the mesons to the quarks are of Yukawa 
type [7(|: 


£qQtt = igQQir1pl 5 T ■ 1Tip, (16) 


FIG. 1: The diagram (a) represents a schematic one-loop dia¬ 
gram of the quark correlator, which can be obtained from the 
two-point function of the viscous-stress tensor for the quark 
constituents. The double-dashed lines for the quark propaga¬ 
tors indicate that they have some finite thermal width which 
can be derived from the quark self-energy diagrams (b) for 
quark-meson loops (for n,a meson). 

Given the meson masses and quark-meson couplings, 
the quark thermal width is obtained from the imaginary 
part of the loop diagrams shown in Fig. [T] (b). The 
off-shell quark self-energy contains four branch cuts on 
the ener gy axis, but on-shell only the Landau cut con¬ 
tributes [691 ] and the result can be written as 

r Q = T Q(Q M ) 

M—tv, a 
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d 3 l 
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jt .,0 , i u ^n Q (l) + n M {u) 
6(k +Wg-W M ) - —j— - 


4 W Q W M 





k°—u 


(23) 


where u = k — l, riM^k) is the Bose-Einstein distribution 


n M {k) 


1 


e^« - 1 


(24) 


with = (fc 2 + m^) 1 / 2 , and 

Lq(q-k) (k, l) = 3 (M 2 - fc • 0 , (25) 

L Q{Qa) {k,l) = ^^(M 2 Q + k-l). (26) 

Note that use of an on-shell quark thermal width is con¬ 
sistent with the quasi-particle approximation we are us¬ 
ing to describe the system. 

The last input needed is the sound velocity c s , which 
is required for the evaluation of the bulk viscosity in 
Eq. (fl5l) . It can be obtained from the pressure P and 
energy density e as: 


c 


2 

s 



(27) 


where s is the entropy density and cy the specific heat 
at constant volume: 


~- r (£) y - <»> 

The pressure P and energy density e are given in the 
quasi-particle approximation to the NJL model as 


P 

e 


= 4N c N f 


d 3 k 

( 2^)3 


n Q (k) 


k 2 /H \ 
/ 


(29) 


We note that the use of the quasi-particle approximation 
for the thermodynamic quantities is consistent with the 
large N c counting [§4j]. Inclusion of meson loops con¬ 
tributions requires care with respect to thermodynamic 
consistency. 


III. NUMERICAL RESULTS AND DISCUSSION 

The free parameters in the NJL model are the cur¬ 
rent quark mass mg, the coupling G and the cutoff 
mass A that is required to regularize vacuum loop di¬ 
agrams. They are fixed to obtain reasonable vacuum 
values for the quark condensate (uu), pion leptonic de¬ 
cay constant / w and and the pion mass m,. With toq = 
m u = rrid = 5 MeV, GA 2 = 2.14 and A = 653 MeV, 
one obtains: (uu) = (—252 MeV) 3 , / w = 94 MeV and 



FIG. 2: (Color on-line) (a) Temperature dependence of Fq 
( solid line) for k = 0.5 GeV and the individual contribu¬ 
tion from ty— (dashed line) and a —quark loops (dotted line), 
(b) Temperature dependence of quark mass (solid line), the 
pion (dotted line), and the sigma (dashed line) meson masses 
and quark-meson couplings (<?qq^/4-7t: dash-dotted line and 
9qq<t/ Jto dash-double-dotted line). The straight vertical red 
line denotes the Mott temperature (Tm = 0.187 GeV). 


= 142 MeV. The vacuum value of the constituent 
quark mass is Mq = M u = Md = 328 MeV mass and of 
the a mass is to ct = 663 MeV. 

We start our numerical discussion from on-shell quark 
thermal width rg(fc,T), which we present in the top 
panel of Fig. (0 for a specific value of the momentum, 
k = 0.5 GeV. Also shown are the individual contributions 
from the 7r— and cr—quark loops. In the lower panel 
of the figure, we show the temperature dependences of 
the quark-meson couplings and quark, it and a masses. 
The straight vertical red line in the figure denotes the 
Mott temperature Tm, at which m^T/vf) = 2 Mq(Tm)- 
For the parameter values used in the present paper, 
Tm = 0.187 GeV. Here one should notice that the nu¬ 
merical strengths of the pion and sigma modes are oppo¬ 
site in sign [36l |70(, but the magnitude of the pion mode 
is approximately three times larger than the magnitude 
of the sigma mode in the temperature range T > Tm- 
In the lower panel of the figure, one should notice that 
at sufficiently high temperatures, the gQQ v ~ gQQa as 
well as m a ~ m w , a feature due to the chiral symmetry 
restoration at high temperatures. 

It is instructive to examine the invariant-mass M de¬ 
pendence of the individual contributions Tq^q^ and 
rQ(Q< 7 ) to Tq; M is given by M 2 = (fc 0 ) 2 - fc 2 . Re¬ 
sults for k = 0.5 GeV and three different tempera¬ 
tures are shown in Fig. [3] The Landau-cut regions 
for quark-pion (Qn) and quark-sigma ( Qa ) loops of the 
quark self-energy are respectively 0 < M < |m w — Mq\ 
and 0 < M < \m a — Mq |; they are identified by 
the nonzero values of the corresponding Fqqm(A7)’s, 
as indicated in the figure. Clearly, the Landau-cut is 
very sensitive to the temperature as it is determined by 
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the temperature-dependent quark mass Mq and meson 
masses m n and m a . The strengths of the Tqqm are 
essentially controlled by the quark-meson coupling con¬ 
stants; please note the different scales in the vertical axes 
of the different panels. Also, while the threshold condi¬ 
tion ( m a — 2Mq > 0) for the Qa loop is satisfied for 
all temperatures, the corresponding threshold condition 
for the Qn loop (m„ — 2 Mq > 0) is only satisfied beyond 
the Mott temperature Tm = 0.187 GeV. This last feature 
explains the fact that the quark-mass pole, indicated by 
the vertical dashed lines in the figure, lies outside of the 
Landau cut for Qn loop at T = 0.180 GeV. 



FIG. 3: Invariant mass distribution of the quark thermal 
width from the Qn (solid line) and Qa (dotted line) loops 
for |fc| = 500 MeV for temperatures (a) T = 0.180 GeV, 
(b) 0.190 GeV, and (c) 0.200 GeV. The curve for r q(q&)(M) 
as been multiplied by — 1 . 

Next we examine the temperature dependence of Tq 
and of it inverse, the collisional time tq = 1/Tq, for three 
different values of momentum; the results are shown in 
Fig. U The peak position and strength of the temper¬ 
ature dependence of Tq is strongly momentum depen¬ 
dent, a feature that reflects the absorption and emission 
processes of the quark interacting with mesonic modes; 
recall that Tqqm physically means the on-shell proba¬ 
bility of forward and inverse scattering between a quark 
Q and the mesonic modes M 0. In forward scatter¬ 
ing, Q may be transformed to a thermalized M by ab¬ 
sorbing a thermalized anti-quark Q, while in the inverse 
process, an off-equilibrated Q may be regenerated via 
M —► QQ dissociation. We should note as well that Tq 
decreases and ultimately goes to zero for large tempera¬ 
tures; the temperature, for which it tends to be numeri¬ 
cally zero, depends of course on k and is determined by 
the temperature-dependent quark and meson masses and 
couplings. The implication of Tq decreasing with T for 
large T has the implication that 77 and £ will increase 
with T, and again, the value of T where it starts to in¬ 
crease and the rate of increase depend on the parameters 
of the model. 

When considering the collisional time, shown in the 



FIG. 4: (Color on-line) (a) Temperature dependence of the 
quark thermal width Tq and (b) collisional time tq = 1 /Tq 
for three different values of momentum: k = 0 (dashed line), 
0.5 GeV (solid line), and 1 GeV (dotted line). The vertical 
red line indicates the Mott temperature. 


lower panel of Fig. |4j one can compare time scales of dis¬ 
sipation with the typical evolution time At of the mat¬ 
ter produced in a HIC, which is At ~ 1 — 10 fm. One 
can clearly identify well-defined temperature ranges for 
which tq < At, i.e. for which the propagating quark 
in the medium suffers sufficient dissipation by scattering 
with the thermalized constituents of the medium. 



FIG. 5: (Color on-line) (a) Temperature dependence of 77 
and (b) rj/s for different values of constant collisional time: 
tq = 5 fm (dashed line), 1 fm (solid line), and 0.2 fm (dot¬ 
ted line). The straight horizontal red line indicates the KSS 
bound, rj/s = 1/An. 

Before using the full temperature- and momentum- 
dependent rQ(fc) in Eqs. ITU) and (fl5l) . we have cal¬ 
culated the temperature dependence of 77 and £ for a 
constant value of tq = 1 /Tq. Although the tempera¬ 
ture and momentum dependence of Fq will modify some¬ 
what the results for the viscosities, the calculation with 
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a constant T q will bring insight regarding the integrands 
in Eqs. © and m when r Q can be taken out of the 
integrals. In Fig. Owe show the corresponding results for 
77 and 77 /s, for different values of tq. While 77 is a mono- 
tonically increasing function of T, the ratio 77 /s exhibits 
two different rates of increase in two different temper¬ 
ature domains, which are associated with two different 
phases. 



FIG. 6: (a) The temperature dependence of entropy density, 
normalized by Stephan-Boltzmann (SB) limiting value ssb, 
given in Eq. (13011 . (b) The squared speed of sound c 2 as a 
function of the temperature. Straight horizontal dotted lines 
stand for the corresponding SB limits. 

The observed change in the rate of increase of 77 /s with 
the temperature can be understood from the temperature 
dependence of s, normalized by the Stephan-Boltzmann 
(SB) limiting value: 

< 30 > 

One observes in Fig. [ 6 ])a) that the rate of increment in 
s/ssb mainly changes around T « 0.175 — 185 GeV; 
above this temperature, the entropy density of quark 
matter becomes almost identical with the SB limiting 
value, which is denoted by straight horizontal dotted line 
in the figure. It is important to note that we are tak¬ 
ing into account the contributions of the quarks only in 
the expression for ssb', including gluon degrees of free¬ 
dom will increase ssb by a factor roughly equal to 1.75 
(for Nf = 2). Being proportional to the slope of s, the 
specific heat cy is magnifying the transition by exhibit¬ 
ing a smooth hump structure around T ps 0.175 — 185 
GeV as shown in Fig. [U])a). Such a behavior of cy has 
been obtained also in other recent calculations employ¬ 
ing the NJL model, as e.g. in Refs. (34, 72l. l73l|. Note, 
however, that the maximum of cy is still well below the 
corresponding SB limit in QCD which, of course, includes 
gluon degrees of freedom and is thereby 1.75 larger than 
the pure-quark SB limit. This is a welcome feature as 
lattice QCD simulations 01 show no indication that the 
specific heat exceeds the QCD SB limit at any tempera¬ 
ture. For calculating the bulk viscosity, we have obtained 


Cg using Eq. (Ell); its temperature dependence is shown 
in Fig. Hb). Assuming the system is not very far from 
equilibrium, all thermodynamical quantities are obtained 
for non-interacting system of quark matter, although a 
finite (not zero) probability of quark-meson interaction 
has to be considered for getting a non-divergent values 
of transport coefficients. 

Fig. [T] shows results for the temperature dependence 
of £ and £/s, for different constant values of tq. A non¬ 
monotonic temperature dependence of £ is obtained, with 
a distinctive two-peak structure at temperatures T « 
0.150 - 0.160 GeV and T « 0.210 - 0.220 GeV. Similar 
double-peak structure for £ was obtained in Ref. [58|. 
Note that the first peak is located at a temperature below 
the Mott temperature and will not be seen (Fig. [TUI) when 
using the full temperature and momentum dependence of 

r Q (fc,T). 



T (GeV) 

FIG. 7: (a) Temperature dependence of / and (b) (/s for 
different values of constant tq . In the lower panel, also shown 
are results from the literature: Ref. I4h| (filled circles and 
squares), Ref. [54| (open circles), Ref. [n| (stars). 

The double-peaked structure is clearly the result of 
a competition between the two conformality break¬ 
ing terms in the integrand of Eq. (fl5l) . 1/3 — Cg and 
d(P 2 Mq) / d/3 2 : the first is associated with Cg ^ 1/3 
and the second is associated with the bulk mass trans¬ 
port (28|. To get a graphical understanding of the compe¬ 
tition, let us define auxiliary quantities containing these 
terms: A = 1/3 — Cg, B = d(f3 2 Mq) /dj3 2 , C = Ak 2 , 
D = Cg B and E = (C — D) 2 . These quantities are plot¬ 
ted in the different panels in Fig. [5] 

The key feature here is that the C (with any fixed 
momentum) and D curves intersect each other at three 
different points; therefore, E = (C — D ) 2 possesses three 
nodes (minima) and, of course, two maxima along the 
T-axis. These two maxima remain when E = E{k) is 
integrated over k in Eq. m and the two-peak structure 
is clearly understood in these terms. Of course, if Mq 
is temperature independent, then one has only one peak, 
as the lower peaks in Figs. [7] and [SJc) will not appear. 






























FIG. 8: Temperature dependence of the auxiliary quantities 
defined in the text. Panel (a): A = 1/3 — c 2 (solid line) 
and B = d([3 2 MQ)/d/3 2 (dashed line). Panel (b): C = Ak 2 , 
for |fc| = 0.5 GeV (solid line) and |fc| = 0.3 GeV (dashed 
line), D = c 2 B (dotted line). Panel (c): E = (C — D) 2 for 
|fc| = 0.5 GeV (solid line) and |fc| = 0.3 GeV (dashed line). 


FIG. 9: (Color on-line) (a) Temperature dependence of r/ and 
(b) rj/s using the full momentum dependence of the quark 
thermal width P g(k). Comparison with results from the 
literature: Ref. [2f| (filled circles), Ref. (4^1 (open circles), 
Ref. [34[ (open squares), Ref. [2f| (filled squares), Ref. [32 | 
(triangles), Ref. [3(j] (stars). 


Likewise, as seen in the lower panel of Fig. [3 the ratio 
C/s also presents a non-monotonic behavior. For the sake 
of comparison, Fig. [7Kb) also displays results from the 
literature. 

Finally, we present results for the viscosities and vis¬ 
cosity to entropy density ratios with the full momentum 
dependence of the quark thermal width TQ(k). It is im¬ 
portant to reiterate that our results have physical sig¬ 
nificance only within a limited range of temperatures, 
approximately Tm < T < 0.240 GeV. The lower limit is 
due to the fact that in the NJL model, when using the 
on-shell thermal width, finite viscosities are obtained for 
temperatures T larger than the Mott temperature Tm 
because Tq(A:) is nonzero only for T > Tm in this case. 
As discussed earlier, the off-shell quark self-energy con¬ 
tains four branch cuts on the energy axis, but on-shell 
only the Landau cut contributes. The upper limit is not 
a precise limit, but it is related to the lack of asymptotic 
freedom in the NJL model; it is known that the large 
temperature behavior of viscosities in QCD, particularly 
of bulk viscosity, is very different from the one derived 
from a non-gauge, contact-interaction model 0 . Given 
this proviso, we discuss next our results for the temper¬ 
ature dependence of 77 and 77 /s, and C, and ( '/s; they are 
presented respectively in Figs. [9] and [10] 

First of all, one notices in Figs. [9] and [TO] that the 
use of the full temperature dependence of Fq(/c) restricts 
the domain of temperatures where finite viscosities are 
obtained. Our numerical values of 77 and 77 /s within the 
temperature range Tm < T < 0.240 GeV, are quite close 
to some earlier estimates in Refs. [HI, [HI, [H, Hj, [H, [45[ . 
Similarly, Fig. m shows fair agreement of our calcula¬ 
tion of C/s with some earlier results of Refs. [H, HU, [ 53 ], 
@. 

The interesting two-peak structure in C shown in 



FIG. 10: (a) Temperature dependence of C and (b) C/s using 
the full momentum dependence of the quark thermal width 
F Q (fc). Comparison with results from the literature: Ref. 53 
(circles), Ref. [2^ (squares), Ref. [53|] (triangles), Ref. 134] 
(stars). 


Fig. m is not visible here because, as discussed above, 
the first peak is below the Mott temperature. In addi¬ 
tion, the figures show that the viscosities have a rapid 
increase with temperature for T > 0.235 GeV. This is 
due to the fact that for high temperatures, the contribu¬ 
tions from quark-meson fluctuations to Tq(k) decrease 
rapidly for large T, as shown in the top panel of Fig. [2] At 
higher temperatures, quark-quark and quark-gluon scat¬ 
terings will become more important than quark-meson 
scattering. To include such processes, one would proba¬ 
bly need a model that incorporates asymptotic freedom, 
as already mentioned in the previous paragraph. As our 
focus of interest in the present is on the contributions 
from quark-meson thermal fluctuations, we reserve for a 
future publication the investigation of the contributions 









































of quark-quark and quark-gluon processes to viscosities. 



T (GeV) 

FIG. 11: r]/s (a) and (/s (b) at high temperature domain. 
HTL (circles) and FRG (stars) estimations of 77 /s by Arnold 
et al. [9] and Christiansen et al. ;37| are compared with 77 /s 
of Eq. (1141) at tq = 4.5 fm (dashed line) and tq = 0.75 fm 
(solid line). While HTL (circles) estimation of r]/s by Arnold 
et al. 0 is compared with C/s of Eq. (1151) at tq = 0.3 fm 
(dashed line) and tq = 5 fm (solid line). 

We conclude with a discussion of the high-temperature 
behavior of 77 /s and C/s by contrasting results from the 
literature and what one would obtain from our Eqs. ED 
and ED when using finite values for the collisional 
time tq. Fig. (fill) is showing results for these ratios 
for T > 0.240 GeV from hard thermal loop (HTL) cal¬ 
culations [9!;, HR Irfl ] and from a functional renormaliza¬ 
tion group (FRG) calculation for pure Yang-Mills the¬ 
ory [37|. For 77 /s we note that the HTL and FRG results 
are markedly different. The FRG result is almost identi¬ 
cal to the one obtained from Eq. ([Till with tq = 0.75 fm 
are seem to match very well with our result shown in 
Fig. ©(b) form temperatures around T = 0.23 GeV. On 
the other hand, for C/s there are qualitative and quanti¬ 
tive difference between the HTL results from those from 
Eq. (fl5l) for the two values of tq used—there are no FRG 
results available for C/s. The discrepancies are significant 
and clearly further studies are required. 

IV. SUMMARY AND PERSPECTIVES 

In summary, we have investigated the temperature de¬ 
pendence of the quark-matter shear 77 and bulk ( viscosi¬ 
ties and their ratios to the entropy density s. The focus 
of interest was on the contributions from Qt and Qa fluc¬ 
tuations to the on-shell quark thermal width Tq, a cru¬ 
cial ingredient in the calculation of viscosities. We calcu¬ 
lated rQ(fc, T) from the imaginary part of the quark self¬ 
energy arising from Qt and Qa loops using temperature- 
dependent quark and meson masses and quark-meson 
couplings. We have investigated in detail the Landau-cut 
structure of the quark self-energy using real-time thermal 


field-theory. The temperature dependence of masses and 
coupling constants were obtained in the NJL model in the 
quasi-particle approximation. The entropy density and 
the speed of sound, which is needed to calculate C, were 
also obtained in the NJL model within the quasi-particle 
approximation. 

The temperature dependence of the masses and quark- 
meson couplings non-trivially influences the Landau-cut 
structure of the quark self-energy. The quark thermal 
width, calculated from the on-shell quark self-energy, be¬ 
comes non-zero and non-negligible within a range of tem¬ 
perature and momentum where the quark pole is situated 
well inside the Landau-cut. When using a temperature- 
and momentum-independent quark thermal width Tq, a 
monotonically increasing shear viscosity as a function of 
temperature was obtained, whereas, the bulk viscosity 
exhibits a two-peak structure near the phase transition 
temperature, reflecting the competition between confor¬ 
mality breaking terms related with the speed of sound 
in the medium and the bulk mass transport. On the 
other hand, when the full temperature- and momentum- 
dependent thermal width is used, finite viscosities are ob¬ 
tained only for temperatures larger than the Mott tem¬ 
perature Tm = 0.187 GeV. In addition, for high tem¬ 
peratures, due to rapidly increasing mesonic masses and 
decreasing of the constituent quark mass, the probability 
for quark-meson fluctuations become negligible and the 
viscosities increase considerable and other processes not 
related to quark-meson processes take place. Between 
the two extremes, our results for the 77 /s and (,/s are in 
fair agreement with results from the literature, although 
marked differences can be observed. 

Future work includes the investigation of viscosities 
at high temperatures by the introduction of quark- 
quark and quark-gluon processes. Of course, agreement 
with standard perturbative QCD results are expected 
when using interactions that incorporate the property of 
asymptotic freedom of QCD. At low temperatures, when 
the medium is dominated by pions and low-mass reso¬ 
nances [38|, the NJL model can be used to obtain pure 
meson-meson interactions, deduced e.g. via the bosoniza- 
tion technique of NJL model, and results are expected to 
be very close to the standard results for the viscosities of 
hadronic matter. 
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